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Article Info ABSTRACT 

Article history: Impairment of cognitive and working memory after stroke was common. 
Vascular dementia (VaD) was a prevalent type of dementia that was caused 

Received Jun 3, 2019 by an impaired blood supply to the brain because of a series of small strokes. 

Revised Aug 6, 2019 Electroencephalogram (EEG) gives information about brain status and 

Accepted Aug 25, 2019 activity, so it had a lot of potential to be used in diagnosing people with 


dementia. Since the EEG signal is extremely non-linear and non-stationary 
data, traditional Fourier analysis such as Fast Fourier Transform (FFT) that 
Keywords: broadens sinusoidal signals cannot describe the amplitude contribution of 
each frequency value in specific time. Meanwhile, Hilbert Huang Transform 
(HHT) was based on the characteristic local time scale of the signal, it can 
efficiently obtain instantaneous frequency and instantaneous amplitude for 
Vascular dementia nonstationary and nonlinear data. In this paper, HHT was employed as 
feature extraction method to extract the energy features of frequency bands 
from post stroke patients and healthy subjects. The extracted features were 
fed into extreme learning machine (ELM) for classifying post stroke patient 
with VaD and healthy subjects. The results of classification accuracy using 
HHT as feature extractor and FFT as feature extractor were compared. 
The mean accuracy of classification using HHT was 59.14%, respectively, 
while mean accuracy of classification using FFT was 94.4%, respectively, 
in classifying post stroke patient with VaD and healthy subjects. 
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1, INTRODUCTION 

Stroke is regarded as the leading risk factor for cerebrovascular disease which affect severe physical 
disability and cognitive impairment [1]. About 30 percent of people with stroke, experience dementia in the 
first year of stroke diagnosis [2]. The precondition for VaD diagnosis is the dysfunction of post-stroke brain 
memory [3]. Following Alzheimer's disease (AD), VaD is the second most frequent type of dementia and 
considered resulting from cerebral vascular disease and ischemic or hemorrhagic brain injury-related [4]. 
Between | and 4% of elderly people aged 65 years tend to suffer vascular dementia and prevalently doubled 
5 to 10 years after this age [5]. 

Due to post stroke, working memory (WM) is considerably impacted in this condition. 
WM provides temporary storage and manipulates information within a short time for complex cognitive 
functions such as attention, comprehension, reasoning, planning and learning [6]. An individual's WM ability 
can momentarily store about 7+2 items [3, 6]. 
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Neuroimaging and non-neuroimaging methods have been used to diagnose stroke. These methods 
play a significant part in the diagnosis of stroke disease. Magnetic Resonance Imaging (MRI), Positron 
Emission Tomography (PET), and Single Photon Emission Computed Tomography (SPECT) are 
neuroimaging methods that widely known used for stroke diagnosis. Although neuroimaging methods are 
commonly used for stroke diagnosis, they still have several disadvantages [7]. The prevalent disadvantages of 
the neuroimaging method include high price, time consuming, and radiation risk imposition. In addition to 
these methods discussed, Electroencephalogram (EEG) is one of the helpful instruments for detecting and 
tracking impacted brain after stroke dementia. This study describes a method of evaluating EEG signals of 
post-stroke patient working memory compared to healthy subjects. 

EEG has been commonly used in clinical diagnosis because of its high temporal resolution and has 
been recognized as an efficient physiological technique that represents the hidden cortical abnormalities by 
offering a measurable understanding for diagnosing or evaluating future dementia severity [5]. 
The abnormalities in the EEG spectrum of dementia patients such as the slowing of EEG signal activity in 
dementia patient can be recognize by using spectral analysis [5]. These changes shown in the power spectrum 
changes in low frequency, in which low frequency power was increased and high frequency power was 
decreased [8]. 

For the study of EEG signals, many signal analysis and processing methods have recently been 
suggested [9-14]. Traditional Fourier spectral analysis was used among these methods to extract EEG signals 
for dementia detection [10]. The Fourier analysis method such as FFT in which broadens sinusoidal signals 
cannot describe the amplitude contribution of each frequency value in specific time [9]. Despite the fact that, 
the EEG is an extremely non-linear and non-stationary process, therefore the Fourier analysis cannot 
precisely represent the frequency part of the signal. To tackle this issue, some time-frequency analysis 
techniques have been suggested and used to obtain EEG features [11]. The short time Fourier transform 
extracts several frames of the signal to be analyzed with a window that moves with time and used to convert 
EEG fragments into EEG spectrogram representations in time-frequency domain [15]. The wavelet analysis 
is comparable to the short-term Fourier Transform which used a variable window size across the entire EEG 
signal length [16]. In fact, the above described techniques of time-frequency analysis have its own advantage 
in non-stationary signal processing, but they are all based on Fourier method and have not given a unique 
definition of frequency mathematically. 

HHT presents nonlinear frequency correlation or time-spatial characteristics of signal analysis, such 
as in ECG and EKG. Besides, HHT is an adaptive and intuitive method which help analyzing the non- 
stationary signals to locate the local changes and providing the ability for time-frequency analysis and the 
adequacy for complex data [15]. Therefore, in this paper, HHT which is designed specifically for analyzing 
non-linear and non-stationary data such EEG signal was used as feature extraction method to extract the 
energy features of frequency bands from post stroke patients and healthy subjects. 


2. RELATED WORKS 

There are different techniques for extracting brain data, such as spatial filtering [17], temporal 
information [18], spectral information and event-related potential (ERP) [19]. Indeed, instead of using 
functional magnetic resonance imaging (fMRI), other researcher used spectral analysis as an alternative 
method [20]. In related study, EEG analysis using spectral-based analysis has been commonly used by late 
researcher until the recent now. Spectral analysis method often been used by other researcher and this method 
provides different opinions in different fields of research and offers useful information. 


2.1. Relative Power Ratio 

Relative Power Ratio (RPR) is used to quantify the EEG signal changes in dementia patients in 
which low frequency is increase in the power spectrum and high frequency is decrease in the power spectrum 
[21]. The RP for each chosen frequency band was computed using (1) [17]. 


¥! Selected frequency range 


0 = 
RP( %) >» Total range (0.5—-64Hz) 


(1) 


2.2. Spectral Entropy 

Previous studies investigated the EEG signal pattern of stroke patients cognitive and thinking 
ability. Spectral Entropy (SpecEn) is used to evaluate signal spectrum flatness and is regarded a helpful way 
to quantify frequency slowdown due to dementia[16]. SpecEn 1s calculated to apply the Shannon's entropy to 
the normalized Power Spectral Density (PSD) as shown in [22]. 
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SpecEn = eater log[PSD,, (f)] (2) 


2.3. Fourier Analysis 
Using spectral analysis is the common method of assessing EEG. Power spectral shows the EEG 


signal frequency content or the frequency power distribution. Power spectrum is really helpful for 
distinguishing activity in the EEG signal [23]. FFT is used to calculate the global EEG value by using Power 
Spectral Density [24]. 


3. METHODOLOGY 
This part addresses this research's experimental procedure from data acquisition to results 


assessment. Figure 1 showing the research's workflow to analyze stroke patient processed data. 
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Figure 1. Block diagram of the proposed framework 


3.1. Inclusion and Exclusion Parameter 

EEG data sets for 35 healthy subjects were recorded, age range 52 years, and 50 post stroke patients 
age ranged 62 years who have severe post stroke after 6 months on set. The patients were recruited from 
Pusat Perubatan Universiti Kebangsaan Malaysia (PPUKM). The ethical standard of the HUKM Committee 
on Institutional and Research Ethics complied with all processes conducted in this study involving human 
respondents. For this study, it took official permission from patients to continue with the EEG data collection 
by filling out the consent form. As shown in Figure 2, EEG data was recorded using BR321 with a total of 32 
electrodes, plus system reference electrodes and the ground with sampling rate of 256Hz and SOHz 
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Figure 2. 32 channels of BRi32 


In this study, 26 channels were used from 32 channels, Fpl, Fp2 AF3, AF4, F7, F8, FT7, FC3, FC4, 
FT8, T7, C3, C4, T8, TP7, CP3, CPZ, CP4, TP8, P7, P3, P4, P8, Ol and O2. The EEG raw data was sampled 
down to 128Hz and filtered from 0.5Hz to 64Hz with bandpass. Then HHT was used as feature extractor to 
extract the feature of EEG data from post stroke patients and healthy subjects. 

Figure 3 shows the time length of data acquisition consists four memorize task and four recall tasks 
for working memory. Each memorize and recall task was conducted for 10 seconds with 5 seconds rest 
between tasks. 
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Total time: 150 seconds — 
Ready time: 5 seconds 


Memorize Task time: 10 seconds 
Recall Task time: 10 seconds 
Rest time: 5 seconds 


End time: 0 - S seconds 





Figure 3. Data acquisition time length experiment 


The subjects are asked to be motionless as possible. Five memorize task and five recall tasks 
appeared on the TV screen alternately for 10 seconds, the subjects were asked to rest for 5 seconds after each 
task. Memorize phase and recall phase was involved in the EEG data collection. EEG signal at recall phase 
will be used as a marker to find out the abnormalities of the EEG signal between healthy person and post 
stroke patient. 


3.2. Hilbert Huang Transform 

In earlier studies, spectral entropy and energy features of Hilbert Marginal Spectrum (HMS) 
correlate to EEG frequency for detection of seizures has been used [25]. 

The HHT is derived from empirical mode decomposition (EMD) and Hilbert spectral analysis [26]. 
EMD 1s decompose the data to get the intrinsic mode function (IMFs) of the signal which are amplitude and 
frequency modulated signal. 

IMF is a data embedded in oscillation mode [26, 27]. Each element must comply with certain 
criteria: the mean value of the two envelopes defined respectively by local maxima and local minima must be 
zero and the maxima and minima value of zero-crossing must be equal or differ to one [26, 27]. Figure 4 
shows how the original signal is decomposed to IMFs. 


Residue = S(t) 

X1(t) = Residue 

valuei=1,k=1 

while Residue =/= 0 or Residue not monotone 

while X has non-negligible local mean 

Upper(t) = spline local maxima X(t) 
Lower(t) = spline local minima X(t) 
Avg(t) = Ya. Upper(t)+Lower‘(t)) 


X(t) = X;(t) — Avg(t) 
value i++ 


end 

IMF (t) = X;(t) 

Value k++ 

Residue = Residue — IMF, 





Figure 4. Pseudocode decompose original signal into IMFs 


The EMD result was a decomposition of [27]: 
S(t) = ¥,UMF,) + residue (3) 


Hilbert Transform is applied to obtain IMFs can be described as: 
1 
H[x(t)] =x*—=y (4) 


Vascular dementia classification based on hilbert huang transform as feature ... (Mohd Ibrahim Shapiai) 


972 O ISSN: 2502-4752 
Lastly, the marginal spectrum h(w), identified by the Hilbert - Huang spectrum, would be described as: 
h(w) = f, H(w, t)dt (5) 


From each frequency value, the marginal spectrum provides a measure of complete amplitude 
(or power) contribution. [26]. 


4. EXPERIMENT RESULT AND DISCUSSION 

The extracted feature of EEG signal based on HHT and FFT from 43 patients of stroke and 27 
subjects of healthy control were fed into ELM to get the classification performance. The parameter setting 
used for the experiment is shown in Table 1. 


Table 1. Parameter Setting of Experiment 


Technique Parameter 
HHT Energy of Hilbert marginal spectrum 
FFT Power spectral density 
ELM Sigmoidal function 


Classification mode 
150 hidden neurons 


Several experiments were performed on the set of hidden neurons used in ELM and the 150 was 
used because it shows higher accuracy. The Table 2 express overall classification performance between HHT 
and FFT. Based on the Table 2, the percentage of mean testing accuracy using HHT as feature extractor was 
59.14% and the percentage mean of mean testing accuracy using FFT was 94.5%. As shown the mean 
percentage of testing accuracy using FFT is higher than using HHT as feature extractor. 


Table 2. Statistic of the Testing Accuracy 


Statistic Testing Accuracy Testing Accuracy 
(HHT) (FFT) 
Mean 0.5914 (59.14%) 0.9445 (94.45%) 
Maximum 0.7619 1.0000 
Minimum 0.5238 0.7143 
Standard Deviation 0.0958 0.1143 


5. CONCLUSION 

In general, based on previous research, HHT method produces better resolution of data than FFT 
with regard to HHT can obtain better data in the frequency domain and time domain than FFT. In this paper, 
HHT was proposed as the feature extractor to classify between the VaD and healthy subjects. However, 
in this study, the result of HHT testing accuracy is lower in which 59.14% than FFT 94.5%. For further 
study, the process of signal decomposing into several IMF and selection of the IMF will be investigated to 
check whether this parameter effect the testing accuracy. 
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